Chiral Corrections and the Axial Charge of the Delta 
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Chiral corrections to the delta axial charge are determined using heavy baryon chiral perturbation 
theory. Knowledge of this axial coupling is necessary to assess virtual-delta contributions to nucleon 
and delta observables. We give isospin relations useful for a lattice determination of the axial 
coupling. Furthermore we detail partially quenched chiral corrections, which are relevant to address 
partial quenching and/or mixed action errors in lattice calculations of the delta axial charge. 
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Introduction. — Numerical simulations of QCD on 
spacctimc lattices provide a first principles method to 
study non-perturbative regime of QCD For the last 
decade, lattice QCD has made dramatic progress due to 
enlarged computing resources, and advances in numeri- 
cal algorithms. Even with considerable progress, lattice 
simulations arc still restricted to unphysically large quark 
masses, and lattice sizes that are not much larger than 
typical hadronic length scales. Fortunately, low-energy 
hadronic physics is dominated by pion interactions which 
can be studied systematically using chiral perturbation 
theory (xPT). Today lattice data in conjunction with 
xPT enable first principles predictions, and in turn, the 
investigation of the low-energy effective theory. 

Properties of baryons can be addressed systematically 
using xPT @, 0|- There are, however, notable compli- 
cations to this proposal. The number of a priori un- 
known low-energy constants for baryons is larger at next- 
to-leading order compared to those of mesons. The chi- 
ral expansion in powers of to^/A 2 , where is the pion 
mass and A x the chiral symmetry breaking scale, is ac- 
companied by an expansion in m^/M^, where Mn is 
the nucleon mass in the chiral limit. Finally the nearby 
delta-resonances can lead to important contributions. 
For a recent review of delta physics, see jij. The ax- 
ial charge of the nucleon, for example, receives impor- 
tant virtual contributions from pion-delta intermediate 
states because the delta-nucleon mass splitting is about 
the same size as the pion mass, and the axial couplings 
Gan, and Gaa are of order one. We focus on the axial 
charge of the delta, Gaa- The value of this parameter is 
largely unknown for two reasons. Firstly and obviously, 
the short mean lifetime of the delta complicates exper- 
imental extraction of this coupling. Secondly the value 
of Gaa has been inferred from one-loop chiral computa- 
tions of various baryon observables. These computations 
are incomplete, however, as they only roughly estimate 
or completely neglect local contributions from unknown 
low-energy constants. There are too many unknowns to 
extract reliable information about Gaa from chiral com- 
putations alone. Lattice QCD simulations can remedy 



this. 

Early work on the delta using lattice QCD centered 
on electromagnetic moments and delta-to-nucleon elec- 
tromagnetic transitions la, J6|. These calculations have 
been refined recently 0, H @| , by including the effects of 
dynamical quarks, and reaching much lower pion masses 
~ 350 MeV. The axial nucleon-to-delta transition has also 
been studied on the lattice for the first time [3, 11 1. 



The delta axial charge, Gaa, can also be determined us- 
ing lattice QCD. A necessarily component for this study 
is the stability of the delta. Above the decay threshold, 
TOjr > A, where A is the delta-nucleon mass splitting in 
the chiral limit, the delta will be a stable particle on the 
lattice. Its static properties can be calculated and the 
effective theory is used, in turn, to extrapolate down to 
the physical point, where the axial matrix element be- 
comes complex valued. Lighter pion masses too can be 
employed. The reason being that the delta decays via p- 
wave pion emission, and the available momentum modes 
are restrictive enough to keep the delta lattice stabilized. 
Lattice study of the delta at lighter pion masses requires 



more care, especially with volume effects [TJ], but could 
better control the chiral expansion of delta properties. 

In this work, we relate various delta matrix elements 
to the axial coupling Gaa- The connected A ++ matrix 
element provides a convenient starting point for lattice 
simulations. We compute the one-loop chiral corrections 
to the delta axial charge, and investigate its pion mass 
dependence. The modulus |Gaa| is shown to be rela- 
tively stable with respect to chiral corrections. Finally 
we perform the partially quenched chiral computation. 

Delta Axial Matrix Elements. — There are vari- 
ous axial current matrix elements in the quartet of A- 
resonanccs. Several of these can be used at zero mo- 
mentum transfer to define the axial charge, Gaa- The 
remaining choices are then completely determined as a 
product of Gaa and isospin Clebsch-Gordan coefficients. 
To arrive at the conventional definition, we define the ax- 
ial charge through the relation 



<A++|J3 5 | A - 



(A-|J M 3 5 |^ 



Gaa-M m , (1) 
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in which appears the axial current J" 5 = Q^^$T a Q , 
where the isospin generators are given by T a = \r a 1 
with r a as Pauli matrices. The factor M.^ encodes the 
spin structure of the forward matrix element, namely 



U u {P)l^l 5 U u {P), where U V (P) is a Rarita- 



Schwingcr spinor. The definition in Eq. |TJ can be easily 
utilized to calculate Gaa on the lattice. Because differ- 
ences of isosinglet matrix elements vanish, we can write 
Eq. U) as 



(A ++ \u 1 ^ 5 u\A + +) c 



GaaM, 



(2) 



-• — az 



where the subscript denotes only connected quark con- 
tractions. 

For completeness, we specify the other matrix elements 
from which the axial charge can be deduced. The isospin 
changing matrix elements are related by the Wigner- 
Eckart theorem. We consider only AI = +1 for ease, 
and find 



(A++|J+|A+) 
<A+|J+|A°) 
(A°|J+ 5 |A-) 



(3) 



where C is proportional to the common reduced matrix 
element. Using isospin, one can relate the isospin transi- 
tions to differences of AI = matrix elements and thus 
leads to 
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Combining all three of these relations shows that C — 
Gaa- One can then utilize the AI = 1 relations in 
Eq. to determine the axial charge, or the matrix ele- 
ment differences in Eq. The former are directly tied 
to weak interaction phenomenology, while the latter are 
straightforward to implement on the lattice, e.g., Eq. $Z§. 
As with the other baryon axial couplings, there are no 
disconnected diagrams to evaluate. When we general- 
ize to partially quenched theories (or additionally mixed 
lattice actions), the above relations between matrix ele- 
ments remain valid because of the vector isospin symme- 
try in the valence sector. 

Chiral Computations.— The SU(2) L ® SU(2) R 
symmetry of two-flavor QCD is spontaneously broken 
down to the vector subgroup. The low-energy dynamics 
are described by pseudo-Goldstone pions emerging from 
spontaneous chiral symmetry breaking. These modes, 



FIG. 1: One-loop diagrams which contribute non-analytic 
terms to the axial charge of the delta. Mesons are repre- 
sented by a dashed line, while the single and double lines are 
the symbols for nucleons and deltas respectively. The solid 
circle is an insertion of the axial current operator. The wave 
function renormalization diagrams are shown at the bottom. 



<j>, arc non-lincarly realized in the coset field E = £ 2 
cxp(2i0//), where 



7T°/V2 



(5) 



and / = 132 MeV is the pion decay constant. Pion dy- 
namics are described at leading order 1 by the effective 
Lagrangian 



f 2 f 2 ?7J 2 

C = ^-tr(0"Et0 M E) + J — ^tr(E f 
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E). 



(6) 



The baryons are contained in SU(2)v multiplcts: a 
doublet N of spin- 1/2 nucleons and a quartet T of spin- 
3/2 deltas. The Lagrangian up to NLO describing the 
nucleons, deltas and their interactions with pions is given 
by 

C = iNv ■ DN — VT^v ■ DT^ + AT^T^ + 2qaNS ■ AN 
+ 9 an (T^N + NA^) + 2g AA T fl S ■ AT", (7) 

where is the baryon velocity, and 5 M the spin oper- 
ator, see 0, Q for further details. The leading order 
axial current derived from Eq. {7} produces the result 
Gaa = 3aa ■ Beyond this order there is a local contribu- 
tion from the NLO current (which only differs from the 
LO current by an insertion of the quark mass) as well as 
loop contributions depicted in Fig. [TJ Evaluating these 



1 Here we adopt the standard power counting: d ~ m, 
e is a small parameter. 
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FIG. 3: Additional one-loop diagrams contributing to the 
delta axial charge in PQxPT. The cross on the loop meson is 
a hairpin contribution. 



FIG. 2: Dependence of the delta axial charge Gaa on the 
pion mass. The plotted bands account for uncertainty arising 
from the unknown NLO coupling. 



contributions, wc find the delta axial charge 
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Ami 



(8) 



The non-analytic functions appearing above, namely, 
C(m,n), J(yn, A, /i), and /C(m, A, \j) are given in [l3|. 
The delta wavefunction renormalization Za appears in 
Lastly the constant A is the parameter appearing 
in the NLO delta axial current. In Fig. [2J we plot the 
pion mass dependence of the axial charge Gaa- Wc fix 
the tree-level axial couplings at their SU(4) values, and 
renormalize the loop graphs so that Rc(Gaa) = 9aa in 
the chiral limit. To qualitatively understand the contri- 
butions from the NLO coupling A, we vary the renor- 
malization scale [i from 800 to 1200 MeV and further plot 
the real and imaginary parts of Gaa along with (minus) 
the modulus. The figure shows that the modulus is gov- 
erned by perturbative chiral corrections, not the real or 
imaginary parts. 

Partially Quenched Chiral Computations. — The 
partially quenched generalization of the Lagrangian in 
Eq. © reads [lj] 



C 



r- 



str 



mQ(S t 



(9) 



where the mass matrix has been generalized to 
niQ = diag(m, m, rrij, rrij, m, m), 



(10) 



in the isospin limit of the valence and sea sectors. The 
valence pion mass m w is given by m 2 = 2Am, while the 
sea pion mass m^- = 2Xrrij . Finally the valence-sea pion 
mass nij u is given by the average of the two. 2 The flavor 
singlet field $o appearing in Eq. (J9]) is $o = strE/v^ and 
has been retained as a device. The axial anomaly allows 
the singlet mass parameter /iq to be on the order of the 
chiral symmetry breaking scale. Subsequently integrat- 
ing out the singlet yields the correct interactions in the 
flavor neutral sector of the theory. These are described 
by the so-called hairpin vertex, see [1^ . 

Baryons fields are described in terms of the 70- 
dimensional supcr-multiplet Bijk containing spin-1/2 
baryons, and the 44-dimensional supermultiplet T^ k con- 
taining the spin-3/2 baryons. For the embedding of 
the familiar nucleons and delta into these supcrmulti- 
plets, as well as their free and interaction Lagrangian, 
see pi, 0, ED, HI- At leading order, the PQxPT delta 
axial current is given by 22 1 



ja 
J fi5 



2H(T S^ + % 



(11) 



where r 



\{£r a tf + £, f T a 0, and r a are partially 
quenched extensions of the isospin generators. Because 



2 These statements may be modified in the case of a mixed lat- 
tice action. It is straighforward to take this into account given 
partially quenched expressions for observables, see [l6l. Il7l [l8| . 
For example, the valence-sea meson mass receives an additive 
renormalization because no symmetry relates the valence and 
sea sectors. 
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the axial charge Gaa can be determined from connected 
quark contractions, Eq. ([2]), we follow [23| and choose 
the upper 2x2 block of f a to be the ordinary isospin 
generators and zeros elsewhere. The constant TL can be 
determined from matching, i.e. TL = <7aa- 



Since we work to next-to-lcading order (NLO) in the 
chiral expansion, we further require contributions from 
the NLO axial current. These involve one insertion of 
the quark mass matrix uiq 



hT^^tS^ stT(M + )+t 4 T^S,T>; k str(rl + M + ) 



(12) 



where the mass operator M + is defined by: M + = \ (^mq^ + £?t£q£). When working to tree-level, there are only 
two independent contributions from the NLO current: one is proportional to the valence pion mass squared and the 
other is proportional to the sea pion mass squared. Compared to SU (2) %PT there is thus one additional parameter 
to be determined. At NLO, there are additionally non-analytic contributions arising from the loop diagrams in Figs. [I] 
and [3] Evaluating these contributions, we find the partially quenched delta axial charge 



Gaa = 9aaZa — 



(4t/) 



1 2g AA C(m JU7 n) + y£ A A 



2 4 

-£(m T7/ u) + -£(m JU ,p) + TZ(r] u , r] u , fj.) 
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9 AN 



:gAAlC(rn ju ,-A, /i) + {g A + g^Jim^, -A, p) - (2g A + gi)J{m ju , -A, p) 



Ami +3171%. (13) 



The partially quenched wavefunction renormalization Z A 



appears in 



14} , and the function lZ(<j), (/>', p) which arises 
from hairpins is given in [l3T | . The constants A and B are 
shorthands for linear combinations of coefficients from 
the NLO partially quenched delta axial current, Eq. (fT2")) . 
With Eq. (fl"3)) . one can extrapolate Gaa in both valence 
and sea quark masses. Finally with trivial modifications, 
the mixed action extrapolation can be performed. 
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